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Abstract 

A quantitative study of the field-induced magnetic ordering in TICUCI3 in terms of a Bose-Einstein condensation (BEC) of 
magnons is presented. It is shown that the hitherto proposed simple BEC scenario is in quantitative and qualitative disagreement 
with experiment. It is further shown that even very small Dzyaloshinsky-Moriya interactions or a staggered g tensor component 
of a certain type can change the BEC picture qualitatively. Such terms lead to a nonzero condensate density for all temperatures 
and a gapped quasiparticle spectrum. Including this type of interaction allows us to obtain good agreement with experimental 
data. 
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TICUCI3 is a spin-dimer system with an excitation gap 
A ~ 0.7 meV in zero magnetic field. The dimers in this 
compound are formed by the S — 1/2 spins of the Cu 2+ 
ions and superexchange interactions are mediated by the 
Cl~ ions. The lowest excited state is a triplet of massive 
bosons (magnons). The density n of the boson mode is di- 
rectly related to the magnetization per site m = gptBU, The 
magnon dispersion has been measured by inelastic neutron 
scattering [1] revealing a considerable dispersion in all spa- 
tial directions. TICUCI3 therefore has to be considered as a 
three-dimensional interacting dimer system. 

A magnetic field H causes a Zeeman splitting of the 
triplet with the lowest mode crossing the singlet state at a 
critical field H c = A/g/iB ~ 5.6 T The ground state for 
H > H c becomes a BEC of this low-energy boson. The 
condensation of magnons corresponds to a phase transition 
into a state with long-range antiferromagnetic order perpen- 
dicular to the applied field. In Ref. [2] it has been shown that 
this transition and the overall shape of the magnetization 
curves as a function of temperature are in rough qualitative 
agreement with the BEC picture. However, we find that the 
simple magnon dispersion A + k 2 /2m used in [2] is only 
justified for temperatures I< IK (see inset of Fig. 1) 
well below the experimental temperature range. Therefore 



* Corresponding author: sirker@phys.unsw.edu.au 
Preprint submitted to SCES '04 



the universal power-law dependence between critical den- 
sity and temperature n c ~ T" with a — 3/2 expected 
within the BEC picture will only be observable at extremely 
low T thus offering a simple explanation for the deviations 
observed in [2], Following Ref. [3] we have therefore first 
calculated the real triplet dispersion using a bond opera- 
tor formalism [4]. The starting point is the strong coupling 
ground state s) where each dimer at site i is in singlet con- 
figuration \i, s). One then expands the Hamiltonian in terms 
of operators t\ a which create local triplet excitations with 
spiral index a. In first order a Hamiltonian bilinear in the 
operators t^' is derived. Only the most important exchange 
pathways between the dimers are considered so that the lim- 
ited exchange parameter set can be determined by a fit to 
the measured dispersion [1]. In agreement with [3] we find 
that terms quartic in the triplet operators derived in next or- 
der yield only small corrections. The most important renor- 
malization of the triplet dispersion is due to the hard-core 
constraint for the triplets, which can be taken into account 
by introducing an infinite on-site potential 
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The corresponding two-particle scattering vertex v(k, ui) 
can be calculated exactly in the dilute limit by a summa- 
tion of ladder diagrams [5]. Doing so we find a considerable 
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renormalization of the exchange parameters compared to the 
mean-field treatment in [3]. For the scattering vertex at the 
band minimum q we find vq — v(q o ,0) ~ 10 meV [4]. 
The density of states p(Q) calculated with the renormalized 
dispersion fi^ is shown in Fig. 1. 
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Fig. 1. Density of states p(Q). The inset shows p(O) near the gap and a 
fit (dashed line) p(S7) ~ y/Q — A yielding an effective mass m e fr ~ 0.2 
meV" 1 . 

For magnetic fields H > H c and temperatures T < A it 
is sufficient to take only the lowest triplet mode into account. 
The Hamiltonian is then given by 

Y S t h+g t k'-g t ''tk' (2) 
h k.k'.q 

where eu = flfe — A and jj,q = g/j,B(H — H c ). In TIQ1CI3 
typical magnon densities in magnetic fields H ~ 6 — 7 T are 
of the order n ~ 10~ 3 [2]. This allows us to calculate the 
thermodynamic properties systematically in an expansion in 
n (gas approximation). Practically it turns out to be sufficient 
to treat the interaction in (2) in the one-loop approximation. 
After diagonalizing the resulting bilinear Hamiltonian we 
can directly calculate the Boson density (magnetization) as 
a function of temperature [4] using the density of states 
from Fig. 1. The result is shown in Fig. 2a. Compared to 
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r^' 3 aT c with r^'^a < 0.1 the one-loop approximation is 
reliable and higher order corrections are small. Therefore 
there must be some physics we missed so far. An interaction 
which we ignored and which is known to be important in 
TlCuCb is spin-phonon coupling [6]. Its main effect will be 
a reduction of the bare magnon bandwidth and therefore an 
effective increase of vo- With vq ~ 25 meV we can indeed 
obtain good agreement with experimental data for T >T C , 
however, the calculated magnetizations at low temperatures 
are still 50% too large [4]. In recent ESR measurements [7] 
a direct singlet-triplet transition has been observed which 
suggests that a perturbation of the form 



n' = 17 ( t qo - 1\ 



(3) 



might be present mixing the singlet and the triplet state. Here 
7 is a small parameter. Such a perturbation is equivalent to a 
staggered magnetic field perpendicular to the applied field. 
It can result from a staggered g tensor or Dzyaloshinsky- 
Moriya interactions which are both induced by spin-orbit 
coupling. With this term included there is no longer a sharp 
phase transition from a condensed into a non-condensed 
phase but instead the number of condensed magnons de- 
creases exponentially with temperature. Furthermore there 
is no longer a gapless mode usually characteristic for the 
condensed phase [4], In TICUCI3 the interaction (3) re- 
quires breaking of the inversion symmetry within a dimer. 
This might be due to small static distortions as discussed in 
Ref. [4]. We therefore expect 7 to be quite small. It turns 
out that it is indeed possible to obtain excellent agreement 
with experiment with 7 = 10~ 3 meV as shown in Fig. 2b. 
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Fig. 2. Experimental M(T) (symbols) for H = 5.9, 6.0, ■ • ■ , 7 T taken 
from Ref. [2]. (a) The theoretical M{T) (solid lines) with A = 0.67 meV 
and vq — 9.8 meV. (b) Interaction (3) included with 7 — 10~ meV, 
A = 0.72 meV and v = 27 meV. 

experiment the theoretically calculated magnetizations are 
too large and the T c 's too high although n c ~ T,? 1 yields a 
good fit with the exponent being close to the observed one. 
Finally, we want to discuss reasons for this shortcoming. 
First, apart from a small temperature interval \T — T c \ < 



